Achieving ground-state polar molecular condensates by chainwise atom-molecule 

adiabatic passage 
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We generalize the idea of chainwise stimulated Raman adiabatic passage (STIRAP) [Kuznetsova 
et al. Phys. Rev. A 78, 021402(R) (2008)] to a photoassociation-based chainwise atom-molecule 
system, with the goal of directly converting two-species atomic Bose-Einstein condensates (BEC) into 
a ground polar molecular BEC. We pay particular attention to the intermediate Raman laser fields, 
a control knob inaccessible to the usual three-level model. We find that an appropriate exploration 
of both the intermediate laser fields and the stability property of the atom-molecule STIRAP can 
greatly reduce the power demand on the photoassociation laser, a key concern for STIRAPs starting 
from free atoms due to the small Franck-Condon factor in the free-bound transition. 

PACS numbers: 03.75.Mn, 05.30.Jp, 32.80.Qk 



I. INTRODUCTION 

A condensate of ground polar molecules with large per- 
manent electric dipoles represents a novel state of mat- 
ter with long-range and anisotropic dipole-dipole inter- 
actions that are highly amenable to the manipulation 
by dc and ac microwave fields 1]. As such, creation of 
such a condensate is expected to be celebrated as another 
milestone that promises to greatly spur activities at the 
forefront of physics research, particularly with respect 
to quantum computing and simulation [2i] and precision 
measurement Q. 

The road to molecular condensation is, however, com- 
plicated by the fact that more degrees of freedom are 
needed to describe molecules than atoms. In particular, 
cooling particles by entropy removal, a direct method 
popular with atoms, has so far proved to be unable to 
lower the temperature of molecules down to the regime of 
quantum degeneracy. Thus, most current experimental 
efforts in both homonuclear ^] and heteronuclear [^-11] 
molecules have all taken a different approach exempli- 
fied by the first experimental realization of ground polar 
RbCs molecules [6|, in which molecules are first coher- 
ently created from ultracold atoms by photoassociation 
(PA) [13] , and are then brought down to the lower energy 
state by a coherent laser field (instead of by spontaneous 
decay (ajIUlIi^)- More recently, by applying a single-step 
stimulated Raman adiabatic passage (STIRAP) JJ,] onto 
the weakly bound Feshbach molecules, groups at JILA 
have successfully created an ultracold dense gas of 
polar '*°K*''Rb molecules. 

In such schemes, there is a relatively large energy dif- 
ference between the initial Feshbach and final ground 
molecular states. The former, being close to the dissoci- 
ation limit, is a highly delocalized state, while the latter 
is a tightly bound state. It is then, in principle, difficult 
to locate a single excited state, capable of a large spatial 
overlap integral [or equivalently a good Franck-Condon 
(FC) factor] with both the initial and final states. The 



desire to overcome this obstacle has led to the idea of 
stepwise STIRAP [11 [isl. and more recently to the idea 
of chainwise STIRAP [16|, both of which are based on 
models where additional intermediate states and Raman 
laser fields are introduced to form a chain of A systems 
[H [ll|. In these STIRAPs, the two lower states within 
each sub-A system are far closer in energy than the initial 
and final states, thereby greatly boosting the chance of 
locating an excited state capable of a large FC transition 
to both lower states. In contrast to a stepwise STIRAP, 
which employs a series of STIRAPs to move molecules 
one step at a time down each lower intermediate state, 
a chainwise STIRAP applies a single STIRAP between 
the initial and final lasers to transfer molecules, while 
uses relatively high intermediate cw lasers to keep all the 
lower intermediate states virtually unoccupied. This lat- 
ter method eliminates the opportunities for molecules to 
inherit decoherence from the unstable lower intermediate 
states, and is clearly an improvement over the former as 
far as the ability to preserve the phase-space density is 
concerned. 

The focus of this paper is on the coupled multilevel 
atomic-molecular condensate systems where the role of 
the initial transition is played by photoassociation. (An 
example is provided in Fig. [1] which will be described 
in detail in the next section.) Our goal is to develop a 
generalized chainwise STIRAP founded on the concept 
of atom-molecule dark state, a coherent population trap- 
ping (CPT) superposition between stable ground species 
[19l . |2(J |. This scheme has several attractive proper- 
ties. First, atoms are directly converted into ground 
molecules. Thus, the loss of atoms typically associated 
with the initial preparation of Feshbach molecules 0, Hlj 
is never an issue here. Second, pulses of longer durations 
can be employed to meet the adiabatic condition; we can 
do so because the atom-molecule dark state is far more 
stable than the molecular dark state, where the initial 
state is highly unstable compared to the ground (atom 
or molecule) states. Finally, the use of intermediate lasers 
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FIG. 1: (Color online) A schematic of a chainwise STIRAP. 



presents us with a new control knob inaccessible to typi- 
cal three-level models. It is the purpose of this paper to 
show that an appropriate exploration of both the inter- 
mediate laser fields and the stability property of the dark 
state can greatly reduce the power demand for the PA 
laser needed for an efficient conversion. This along with 
other efforts involving the use of Feshbach resonance [l^l 
may help to combat the weakness in photoassociation, a 
key concern to STIRAPs starting from free atoms due to 
the free-bound FC factor being typically very small. 

Our paper is organized as follows. In Sec. II, we de- 
scribe our model and the underlying mean-field equations 
and provide the rationale that justifies the implementa- 
tion of chainwise STIRAP in our model. In Sec. Ill, 
we derive a set of linearized equations around the time- 
dependent CPT solution and obtain an adiabatic condi- 
tion via a state expansion over a set of orthonormal base 
vectors provided in Appendix B. In Sec. IV, we apply 
this adiabatic theorem to facilitate our numerical studies 
of the examples that serve to illustrate the main physics 
outlined in the previous paragraph. Finally, a summary 
is given in Sec. V. 



II. MODEL AND CPT STATE 

Figure [T] is the energy schematic diagram for a mini- 
mum model that is capable of illustrating all the main 
points that we want to convey in this paper. A laser 
field associates atoms from two distinct species of states 
|0i) and IO2) into molecules of state |1) in the excited 
electronic manifold with a coupling strength ili propor- 
tional to the laser field and the free-bound FC factor. Si- 
multaneously, a series of laser fields of (molecular) Rabi 
frequency fi^ {i > 2) is applied to move the molecules 
from the excited to the ground state |6) via additional 
intermediate energy states. In our notation, a molecular 
state \i) (i = 1, 2, • • • ,6) is coupled to the atomic states 
via an i-photon process characterized with an i-photon 
detuning defined, respectively, as Ai = tui — Ei/h, 

A2 = (Wl - UJ2) - E2/h, A3 = (Wi - U}2+ W3) - Es/h, 

etc., where uji stands for the (temporal) frequency of the 
laser field with Rabi frequency fi^, and Ei for the energy 
of molecular state \i) relative to the free atomic energy 



level. Further, intermediate states \i) (i = 1,2,- •• ,5) 
are assumed to be unstable; for each intermediate state 
a decay rate 7^ is introduced to describe phenomeno- 
logically the loss of its molecules due to various incoher- 
ent processes. 

As a proof of principle, we consider, in this paper, a 
uniform condensate system with a total atom number 
density n and describe such a system with a set of field 
operators ^'i, where is the operator for annihilating a 
bosonic particle in condensate state \i) . By following the 
mean-field treatment of photoassociation at zero temper- 
ature [2^ in which each is treated as a c number ^i, 
we obtain, from the Heisenberg's equations for operators 
^i, a set of coupled Gross-Pitaevskii's equations for the 
normalized condensate fields "0^ — "^i/ ^/n: 

«V'Oi = -^-052 V"!, (la) 

^■^Os = -^"001 "01, (lb) 



«0i = (Ai - 171) ?Ai + ^-001-002 + -^^2, 



i0i = (Aj 



fie 



-0, 



i+l, 



106 = A606 -I- 

where the molecular Rabi frequencies 



Qi = V2y/^n[''^^ (vi\v 



.,(0) 



(Ic) 

2,3,---5,, 
(Id) 

(le) 



(2a) 

= {v^\v,+i) , i = 2,3,---6, (2b) 

are expressed in terms of the mean electronic Rabi fre- 
quency Cl'f^\ the free-bound FC factor and the 
bound-bound FC factor {vi\vi+i), where Vi and v^^"^ are 
the stationary wave functions (of interatomic distance 
R) for a bound molecular state \i) and a pair of atoms in 
states |0i) and IO2), respectively [13, HI]. In arriving at 
Eqs. H]), without the loss of the main physics, we have 
followed Refs. 0, [13 and ignored all the two-body s- 
wave collisions. Further, in order to better illustrate the 
essential physics, we will limit our study to a model in 
which = fis = and 0,2 = = ^e, where subscript 
o and e stand for the intermediate lasers of odd and even 
indices, respectively. 

Before moving ahead, we note that Jaksch et. al. 
[14l] have identified a set of rovibrational levels from 
X^E+ (ground) and (excited) electronic manifolds 

to implement the honionuclear version of the model in 
Fig. [T] for producing Rb2 molecules in the ground state 
X'^Y.+ {v = 0,j = 0) [=|6)]. It is true that in order to 
arrive at a similar set of pathways for the heteronuclear 
model, one must perform a careful analysis of experi- 
mental spectroscopic data and possibly, ab initia calcu- 
lation of various overlap integrals [FC factors defined in 
Eqs. (O] [26 28]. However, the selection rules for het- 
eronuclear molecules are actually more relaxed; the het- 
eronuclear molecular orbitals do not have g/u symmetry. 
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and the number of possible transitions between the free 
atomic and the ground molecular state is thus doubled. 
As a result, it is not difficult to see that such a model can 
be easily generalized from homonuclear to heteronulcear 
molecules. 

By ignoring the decays (see the justification that fol- 
lows) and subjecting our system to the conservation of 
total particle number, 

6 

1=1 

and that of atomic population difference: Voi ~ "002 = 
(or ■0Oi = "002 = 00 fo'^ ^ balanced model), we find 
that under the conditions of two-, four-, and six-photon 
resonance, namely, A2 — A4 — Ag — 0, the system at 
steady state supports a superposition involving all the 
lower states with the following amplitude distribution 
(see Appendix A for a detailed derivation): 

V'o = , (4a) 
yi + 7l + 2(a2^i/r!6)' 

f>2 = -ivl, = fa<Po> V6 = -a^^ifil, (4b) 

"6 

y'l = V'a = <^5 0, (4c) 

where a = r^o/f^e, £, — ^i/^e- In arriving at Eq. Q, we 
have only retained the leading order term in ^, assuming 
that the intermediate fields (fio, ^e) are far stronger than 
the initial and final fields (fii, Jig) (Unless stated 

otherwise, a similar perturbative interpretation applies to 
all the other results.) It needs to be stressed that Eq. @ 
is derived when all the decays are ignored. As a result. 



it represents a steady state (or CPT state) of Eqs. ([T]) 
(where all the decay rates are included) only when it does 
not involve any unstable states. The only unstable states 
in Eq. (U) are the lower intermediate states: |2) and |4) 
whose amplitudes scale as ^. Thus, it is in the limit 
^ <C 1 when states |2) and |4) remain virtually empty 
that this superposition can be truly called a "CPT" or 
"dark" state. Also evident is that in the same limit and 
for a fixed a, the initial and final populations are solely 
determined by the ratio fii/rig- This lays the foundation 
for converting all the atoms into the ground molecules 
by a chainwise STIRAP where a counterintuitive pulse 
sequence is maintained only between the initial and final 
laser fields. Another important feature is that the initial 
and final populations are function of a^Jli so that the 
change by ili can also be accomplished by varying a. 
However, the full impact of a on the STIRAP has to 
wait until we know the adiabatic condition. 



III. ADIABATIC CONDITION 

To obtain the adiabatic condition, we linearize Eqs. 
^ around the instantaneous CPT state [t) accord- 
ing to i/ii {t) — ifi (t) + 5%pi (i), where (5-0^ (t) is the small 
perturbation, and ipi (t) are Eqs. ^ when fli are re- 
placed with their instantaneous values fli (t) at time t. 
This procedure results in a matrix equation for the ket 

\H) = (<^V'0i,i5V'02>'^V'i>--- ,<5-06)^: 

4 - (M - h) \Sib) (5) 
at 

where 
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with n'l — ipo^i- Following the standard procedure 
[29| . we expand ISiIj {t)) = ^ q {t) \Xj {t)) in the space 
spanned by the instantaneous collective modes \Xj (t)) 
defined by 

M\X,{t))=X,{t)\X,it)). (7) 
In this new basis and under the condition that 



(A,| M |Aj) / |Aj - A,;| << 1 (A,; ^ Xj), Eq. ^ becomes 

i-j^ = XiCi - (Ail 7 |Aj) Cj - i {Xi\ip) (8) 
j 

and can be put in a form convenient for us to estimate 
the magnitudes of various Ci due to the time variation 
of the CPT state \(p), or in another words, to arrive at 
the adiabatic condition. (Note that whenever confusion 
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is unlikely, we omit argument t to time-dependent vari- 
ables.) 

Thus, we see that the most crucial step in develop- 
ing an adiabatic theorem is to determine from Eq. ([7]) 
a set of base vectors |Aj) upon which we can expand 
\6ip (t)). As this step itself can often be quite involved, 
we focus on a simplified model with Ai = A3 = A5 = 0. 
As we show in Appendix B, for this special model and 
in the limit ^ ^ 1, we can apply perturbation theory 
to obtain simple analytical solutions from Eq. ([T]). In 
what follows, we simply quote the relevant results and 
refer interested readers to Appendix B for details. In 
a nutshell, the collective modes are found to consist of 
Ai,2 = [same as Ai^2 in Eq. (jB9[) ] with a double degen- 
eracy, a set of "soft" modes = ztfie/// (2/3) [same 
as A3^4 in Eq. (|B9[) ] that scale as {rt[ and Hq), where 

a*, and two 



sets of "stiff' modes As^g 
as A 



and /3 = Vl + 



= ±0.5f^e\/r 



same 



(0) 



A5 6 in Eq. (|B4bj) ] and A7,8 = ±0.5f^eVl - a + a'^ 
[same as Ay^g in Eq. (|B4c|) ] that scale as fte and are 
larger than the "soft" modes by a factor of Each set 
of nondegenerate modes is symmetrically displaced with 
respect to the modes, and the appearance of the neg- 
ative modes is expected because the CPT state is not a 
thermodynamical ground state. 

The mode b eing a doublet distinguishes the het- 
eronuclear model [30| from its homonuclear counterpart 
which only supports a single mode [2^. In general, due 
to the ground intermediate states being unstable, the two 
modes here are both unstable, in contrast to the dark 
state in a typical three-level system, which is a stable 
superposition, completely isolated from other unstable 
states. However, by choosing the two modes in the or- 
thonormal form 



|Ai) 
IA2) 



-i= (-1,1, 0,0, 0,0, 0,0)^, (9a) 
(-rifi, -ne, 0, 2Cne,0, -2<r!6, O, 2a^n[f 



(9b) 



where |Ai) is exact while IA2) is correct up to the first 
order in ^' = fl'^/ile which scales as ^ <C 1 [see the 
discussion bellow Eqs. dU] because ^' = ipo^j we find 
that |Ai) is completely decoupled from other states as in 
a true dark state while IA2) couples to other states with 
strengths that are at the same order of magnitude as its 
decay rate, 

(A2I 7 IA2) = 2 (72 + a'74) ^'^nl/nlff = l/ro. (10) 



The result in Eq. (fTO| corroborates our intuition that the 
use of relatively high intermediate laser fields can indeed 
make the lifetime of our CPT state, tq, far longer than 
those of the lower intermediate states. Clearly, in order 
to justify the use of the CPT state in Eqs. ^ as the 
adiabatic state for a STIRAP process, we must design 



the STIRAP in such a fashion that it is slow compared 
with the periods of the nonzero modes but fast compared 
with To, the lifetime of the CPT state. 

As a result, we consider all the coupling coefficients 
(A2 1 7 I Ai) involving | A2) weak and ignore them from Eqs. 
([8|). In addition, we also ignore all the stiff modes, be- 
cause they are times more difficult to populate than 
the soft modes, whose eigenvectors are given by 

|A3,4) = — (7^,7^,± — ,0,T-,0,±- 



(11) 

Under these conditions, Eqs. ([8|) are simplified into 
(iA3 -f (A3I 7 IA3)) C3 + (A3I 7 IA4) C4 = - (A3|(^) , (12a) 

(A3 I 7 IA4) C3 + {iM + (A4I 7 IA4)) C4 - - (A4|(^) . (12b) 

where we have ignored 63^4 assuming that they are suffi- 
ciently small compared to A3. 4 in the adiabatic limit. 

Finally, with the help of Eqs. ([6]) and ([Tl]) . we find from 
Eqs. p2|l that to a good approximation, the adiabatic 
parameter ^[ 



r=VW+^?/2 (13) 
can be estimated according to 



V2A2 (1 + 2a4^2)i/4 1 ^ ^1 + 2a^x 



« 1, (14) 



where x = ^il^Q-, A = ile/// (2^) , and 7e// 
(aV+a'73 + 75)/(2/32). 



IV. DISCUSSION 

In what follows, we seek to gain from the r value in 
Eq. ([H]) insights into the parameters, especially, a, that 
optimize the final conversion efficiency. In a dynami- 
cal process where the r value varies with time, we find 
the r value evaluated at time ts to be a good figure of 
merit that distinguishes different STIRAPs, where tg is 
the time when about 50% atoms would be converted into 
molecules if the system were to follow the CPT state. 
With ts defined above and the Gaussian pulses defined 
below, we find 



t. = 



J.2 



2 {h - h) 



In 



2^1 



where T, ti^, and fi? 6 9''"^' respectively, the width, 
peak times, and peak strengths of the Gaussian pulses: 
f^i,6(t) = Vl\j^cyi]y[- [t - ti^ef /T'^]. In all the calcu- 
lations, 71,3,5 — 10*" 
te = LOT,' fl? 



s \ 72,4 = 10 s" 



u = 2.5r, 



3.3 X 10^ s-\ 



2 X 10** 



and 



fig = 8 X 10'' s ^ . Note that although in theory high adia- 
baticity can always be gained at the expense of a large fij. 



5 



10 



^ 10 



10 ■ 



'\ 

T=10^s / \ 1 \ 

T=55ns ^"^ ■■' 

- - T=300ns | - 
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d>1 






FIG. 2: (Color online) (a) The r value at ts as a function of a 
under different T . The population dynamics with (b) T = 10 
\is and (c) T = 300 /is for different q. Other parameters are 
defined in the text. 



in practice, rij is quite limited due to the relative weak- 
ness in photoassociation. For this reason, we have chosen 
fi? to be much weaker than ^% w At such 



57]* J find, using n = 10 



20, 



and /i 



X 10 



131, 



(reduced mass of K and Rh atoms), that il^/wp 
where Wp = hn^/^ /2n In a STIRAP process, due to 
quantum interference, the molecular population in state 
|1) remains extremely small so that even when ^l^/tOp is 
in the order of 10^, rogue photodissociation of molecules 
in state |1) is shown to produce a negligible fraction of 
noncondensate atom pairs (32| . As a result, we ignore 
the rogue photodissociation in this work. 

Figure illustrates how the r value at ts changes 
with a under different T. The most interesting feature 
here is that, for a given T, the r value is quite high on 
the side of a < 1 (the left side of the dashed vertical line 
a = 1), but significantly smaller within certain region on 
the side of a > 1. As a result, we see in Fig. [2jb) with 
T = 10 /is that when a changes from 0.5 to 3.8 (at which, 
r value is near minimum), the conversion efficiency rj[= 
2ipQ (oo)] increases from 1.8% to 22.6%, a trend consistent 
with Fig. [Ha) with T = 10 /is. But, we caution that no 
efficiencies significantly higher than 22.6% are possible 
in this case by further raising the a value because the 
r value actually increases with a when a is sufficiently 
large according to Fig. [2ja); we trace this to the fact 
that unlike the CPT state in Eq. A and jeff and 
hence the adiabatic condition in Eq. (|14p do not scale as 



An important point to make is that were the lifetime 
of the dark state limited to the order of 10 iis, a higher 
efficiency would indeed have to come at the expense of 
a higher PA laser power. This, however, is not needed 
owing to another important virtue of our system - the 
stability of our CPT state, whose lifetime can be made 
much longer than those of the lower states in the ground 
electronic manifold. This, therefore, affords us with a 
plenty of room to increase efficiency by using pulses with 
longer durations rather than higher powers. Indeed, the 
population dynamics in Fig. HJc) with a much longer 
pulse {T — 300 lis) demonstrates, in addition to a trend 
same as in Fig. [Hb), a dramatic increase in the maximum 
efficiency, which can now reach more than 85%. 

To estimate the required power on the photoassocia- 
tion laser field, we consider heteronuclear molecules in- 
volving Rb atoms, for example, KRb [III and RbCs 0]. 
The free-bound FC factor for heteronuclear molecules is 
expected to be smaller than that for their homonuclear 
counterparts because the excited potential at large in- 
ternuclear distance R for the heteronuclear molecules is 
dominated by the van der Waals potential (R~^) and 
thus has a shorter range than that for the homonuclear 
molecules, which is dominated by the resonant dipole- 
dipole interaction (R"'^). An encouraging news accord- 
ing to ab initia calculations in Refs. [27l Issj is that 
the former is only slightly smaller than the latter. As 
a result, in our estimation, we choose a free-bound FC 
factor {^wilii^^^) = 3 X 10"^''^ m'^/^ several times smaller 
than that of Rb2 molecules, which can be on the order of 
lO'^^m'^/^ (for n = 10^° m^'^) according to Naidon and 
Masnou-Seeuws ^2§\. Finally, using 1.6 eao as the atomic 
dipole moment [34] with e the electron charge and ag the 
Bohr radius, we estimate the peak PA laser intensity to 
be 3.8 X 10^ W/cm^. This admittedly high (and yet 
attainable) intensity can be put in perspective by com- 
parison with the case where T = 10 /iS and a = 1, where 
an intensity of more than 100 times higher is needed to 
achieve the same level of high efficiency [in Fig. [2][c) with 
a = 3.81. 



V. SUMMARY 

In summary, we have generalized the chainwise STI- 
RAP from pure molecular to coupled atomic-molecular 
systems. In addition to the known advantages, for ex- 
ample, the increased chance to locate pairs of Raman 
transitions with large FC factors, we have uncovered ad- 
ditional virtues. In particular, a, a ratio between interme- 
diate laser fields, was found to serve as a robust experi- 
mental control knob, inaccessible to the usual three-level 
systems. This control knob together with the stability 
of the atom-molecule dark state may bring us one step 
closer to overcome the PA weakness, so that the ground 
polar molecules can be created directly from degenerate 
atomic gases in a manner that preserves the phase-space 
density. 
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Appendix A 

This appendix provides the steps that we take to ar- 
rive at Eqs. dH). We begin with Eqs. ^ at steady state, 
where the derivatives on the left-hand sides are all set to 
zero. Next, we ignore all the decays as well as all the ex- 
cited populations (-01 = V'3 = '05 = 0). We then see that 
the equations for ■02, V'4j ^-nd ^pe lead to the CPT condi- 
tion: A2 — A4 = Ag = 0, while the rest of equations are 
simplified to 

-C^OiV'Os = V'2, (Ala) 

— Q02 — 04, (Alb) 
-^04 = 06. (Ale) 



where as in the main text we have made the use of ft2 = 
rii — rie, = = f2o, ct = ^o/^e and ^ = Qi/rie. 
For a balanced system with 0g^ — 0g^ = 0g, we find 

from Eqs. (jAl[) that 02 = — C0o, '04 = ckC'0oi ^'^d 06 = 
— (f2o/r2g) a^0Q, which, when combined with the particle 
number conservation in Eq. ([3]), gives rise to 



1 + ^1 + 2 (1 + a2) ^2 + 2 {a^rii/ner 
V2 = -^(fl, <P4 = ^aifl, ipe = -a'^^^l- 

"6 

Clearly, we see that up to the first order in ^, Eqs. (jAl[) 
become Eqs. ([4]) in the main text. 



Appendix B 

In this appendix, we show how to obtain from Eq. ([7]) 
the eigenvalues and eigenvectors (that are needed to de- 
rive the adiabatic condition) in the limit of f ^ 1 for the 
special case of Ai,3_5 — 0. To begin, we divide M in Eq. 
© into two parts: 

M = Mo + M' (Bl) 

where 
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Here, Mg is an unperturbed part including all the inter- 
mediate Rabi frequencies fig and Jig, while M', consisting 
of only the initial and final fields fl'^ and rig, can be re- 
garded as a perturbation to Mg in the limit of ^ ^ 1. 
Next, we determine from the equation 



Mg 



l(0) 



(B3) 



'^1,2,3,4 — 



A 



(0) 

5,6 



(g) _ ifleVl - a + a'^ 
^^7,8 - 5 



(B4a) 
(B4b) 
(B4c) 



eigenvalues A^^'' and eigenstates 



(0) 



of the unper- 



turbed part Mg. We find that A^*'"' take the following 
values: 



The zero eigenvalue has a four-fold degeneracy, and the 
corresponding (orthonormalized) eigenstates are found 
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from Eq. (jB3|) to take the form 



the following equation 



r) 



Af) 



(1,0,0,0,0,0,0,0)^, 
= (0,1,0,0,0,0,0,0)^, 
(0,0,0,0,0,0,0,1)^, 



a 



a 



1 



(o,o,:^,o,-^,o,-,o)-. 



(B5a) 
(B5b) 
(B5c) 

(B5d) 



where P = vl + a^ + a^. This degeneracy, however, 
can be partially lifted by the perturbation M' as we will 



show shortly. Ag^'g and Ay"g are proportional to fl^ and 
are far larger in magnitude than those split from the zero 
eigenvalue by M'; the modes associated with the former 
eigenvalues are far more difhcult to populate than those 
associated with the latter eigenvalues. Thus, it sufHces, 
for our purpose, that we only focus on the eigensubspace 
spanned by the four base vectors in Eqs. (IBSp . in which 
M' in Eq. (jB2p has the following matrix representation 



AO) 



M' 



1 



n'^a^ 
n'^a^ 
f^e 



(B6) 



where the use of M' . = ( A 



,(0) 



M' 



Aj°''\ has been made. 



In the spirit of degenerate perturbation theory [35j, 
we form the following linear combination of degenerate 
states in the four-dimensional Hilbert space 



|A„) = 



4 



>(0) 



(B7) 



where a„ = (a„i, a„2, Ons, clua)'^ is the eigenvector of ma- 
trix M' with an eigenvalue A„ or equivalently it satisfies 



M'a„ = A„a„ 



(B8) 



By solving Eq. (|B8[) . we find the following set of eigen- 
values 



Al,2 = 0, A3 4 = ± 



eff 



2/3 



(B9) 



where ^eff = \/2fl'^a^ + fig. As can be seen, M' re- 
duces the degeneracy of the zero mode from four folds to 
two folds, creating a pair of so-called soft modes, whose 
eigenvalues, A3_4, are symmetrically displaced from zero 
eigenvalue. The (unit normalized) eigenvectors a3^4 for 
A3_4 are found from Eq. (|B8|) to take the form 



(BIO) 



which, when combined with Eq. (jB7p . yields the desired 
soft modes in Eq. (fTT|) . 

At this point, we stress that the zero mode of two- fold 
degeneracy cannot be lifted by M' as one can easily check 
directly from Eq. (O with Xj = that it supports two 
linearly independent (but nonorthogonal) solutions 



|6i) = (-1,1,0,0,0,0,0,0)^, 
\b2) = f-^,0,0,l,0,-a,0,-^ 



(Bll) 

(B12) 



Finally, we apply Gram-Schmidt orthogonalization to 
transform |6i,2) into a set of orthonormalized vectors 
|Ai,2) in Eq. (Id. 
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